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Anti-Jamming  Techniques  for  GPS  Receivers 


Executive  Summary 


This  is  the  final  report  for  the  Air  Force  Research  Lab  (AFRL)  contract  no. 
F30602-99-2-0504.  It  provides  research  results  obtained  over  the  period  of  April  1999  to 
January  2001. 

The  contributors  to  the  research  over  the  span  of  the  contract  are  Professor 
Moeness  Amin  (PI),  Dr.  Yimin  Zhang  (Postdoctoral  Fellow),  and  Mr.  Liang  Zhao 
(Graduate  Student)  from  Villanova  University.  Dr.  Alan  Lindsey  (Project  Manager)  from 
AFRL  has  worked  closely  with  the  research  team  at  Villanova  University  and  has 
provided  valuable  insights  into  several  issues  vital  to  progress  and  advances  in  research. 
Both  Dr.  Zhang  and  Mr.  Zhao  continue  to  work  on  GPS  signal  processing  for  ant¬ 
jamming.  They,  along  with  Prof.  Amin,  are  supported  by  a  new  contract  from  the  AFRL. 

The  fundamental  objective  of  this  research  project  is  to  improve  the  GPS  receiver 
performance  when  subjected  to  strong  nonstationary  interference.  The  techniques  applied 
and  proposed  to  achieve  this  objective  rely  on  the  estimation  of  both  time-frequency  and 
spatial  signatures  of  the  jammers. 

The  report  consists  of  four  chapters  addressing  important  problems  in 
nonstationary  interference  mitigation  in  GPS  receivers.  Each  chapter  has  its  own  abstract, 
introduction,  equation  numbers,  figure  numbers  and  captions,  appendices,  conclusions, 
and  references.  The  first  two  chapters  deal  with  a  single  antenna  receiver,  whereas  the 
last  two  chapters  consider  the  presence  of  multi-antenna  array.  In  all  four  chapters,  the 
nonstationary  interference  is  cast  as  an  FM  signal  and  it  is  mitigated  using  its  temporal 
and  spatial  characteristics  through  subspace  projection  methods.  The  first  three  chapters 
are  GPS  specific  and  utilize  the  GPS  signal  structure  and  its  deterministic  nature.  Chapter 
4  applies  to  the  general  problem  of  suppressing  instantaneous  narrowband  signals  in 
broadband  communication  platforms. 

Subspace  projection  techniques  are  applied  in  Chapter  1  as  a  pre-correlation 
signal  processing  method  for  FM  interference  suppressions  in  GPS  receivers.  The  FM 
jammers  are  instantaneous  narrowband  and  have  clear  time-frequency  (t-f)  signatures  that 
are  distinct  from  the  GPS  C/A  spread  spectrum  code.  In  the  proposed  technique,  the 
instantaneous  frequency  (IF)  of  the  jammer  is  estimated  and  used  to  construct  a  rotated 
signal  space  in  which  the  jammer  occupies  one  dimension.  The  anti-jamming  system  is 
implemented  by  projecting  the  received  data  sequence  onto  the  jammer-free  subspace. 
Chapter  1  focuses  on  the  characteristics  of  the  GPS  C/A  code  and  derives  the  signal  to 
interference  and  noise  ratio  (SINR)  of  the  GPS  receiver,  implementing  the  subspace 
projection  techniques. 

Frequency  modulated  signals  in  the  frequency  band  1.217-1.238  GHz  and  1.565- 
1.586  GHz  present  a  source  of  interference  to  the  GPS,  which  should  be  properly 
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mitigated.  The  problem  of  mitigating  periodic  interferers  in  GPS  receivers  using 
subspace  projection  techniques  is  addressed  in  chapter  2.  In  this  Chapter,  the  signal-to- 
interference-and-noise  ratio  (SINR)  of  the  GPS  receiver  implementing  subspace 
projection  techniques  for  suppression  of  FM  jammers  is  derived.  The  general  case  in 
which  the  jammer  may  have  equal  or  different  cycles  than  the  coarse  acquisition  (C/A) 
code  of  the  GPS  signals  is  considered.  It  is  shown  that  the  weak  correlations  between  the 
FM  interference  and  the  gold  codes  allow  effective  interference  cancellation  without 
significant  loss  of  the  desired  signal 

Subspace  array  processing  for  the  suppression  of  FM  jammers  in  GPS  receivers  is 
introduced  in  Chapter  3.  In  this  Chapter,  subspace  projection  array  processing  techniques 
cire  applied  for  suppression  of  frequency  modulated  (FM)  jammers  in  GPS  receivers.  In 
the  proposed  technique,  the  instantaneous  frequency  (IF)  of  the  jammer  is  estimated  and 
used  to  construct  the  jammer  subspace.  With  a  multi-sensor  receiver,  both  spatial  and 
time-frequency  signatures  of  signal  arrivals  are  used  for  effective  interference 
suppression.  Chapter  3  considers  the  deterministic  nature  of  the  GPS  C/A  code.  The 
receiver  SINR  is  derived  and  shown  to  offer  improved  performance  in  strong  interference 
environments. 

Combined  spatial  and  time-frequency  signatures  of  signal  arrivals  at  a  multi¬ 
sensor  array  are  used  in  Chapter  4  for  nonstationary  interference  suppression  in  direct- 
sequence  spread-spectrum  (DS/SS)  communications.  With  random  PN  spreading  code 
and  deterministic  nonstationary  interferers,  the  use  of  antenna  arrays  offers  increased 
DS/SS  signal  dimensionality  relative  to  the  interferers.  Interference  mitigation  through 
spatio-temporal  subspace  projection  technique  leads  to  reduced  DS/SS  signal  distortion 
and  improved  performance  over  the  case  of  a  single  antenna  receiver.  The  angular 
separation  between  the  interference  and  desired  signals  is  shown  to  play  a  fundamental 
role  in  trading  off  the  contribution  of  the  spatial  and  time-frequency  signatures  to  the 
interference  mitigation  process.  The  expressions  of  the  receiver  SINR  implementing 
subspace  projections  are  derived  and  numerical  results  are  provided. 

Following  this  executive  summary  are  the  lists  of  papers  submitted/published 
based  on  the  above  contributions.  The  four  chapters  included  in  this  report  properly 
integrate  and  extensively  cover  the  material  presented  in  all  paper  submissions  under  this 
contract. 
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Chapter  1 


Subspace  Projection  Techniques  for  Anti-FM  Jamming 

GPS  Receivers 


Abstract 

This  report  applies  subspace  projection  techniques  as  a  pre-correlation  signal  processing 
method  for  the  FM  interference  suppressions  in  GPS  receivers.  The  FM  jammers  are 
instantaneous  narrowband  and  have  clear  time-frequency  (t-f)  signatures  that  are  distinct 
from  the  GPS  C/A  spread  spectrum  code.  In  the  proposed  technique,  the  instantaneous 
frequency  (IF)  of  the  jammer  is  estimated  and  used  to  construct  a  rotated  signal  space  in 
which  the  jammer  occupies  one  dimension.  The  anti-jamming  system  is  implemented  by 
projecting  the  received  sequence  onto  the  jammer-free  subspace.  This  report  focuses  on 
the  characteristics  of  the  GPS  C/A  code  and  derives  the  signal  to  interference  and  noise 
ratio  (SINR)  of  the  GPS  receivers  implementing  the  subspace  projection  techniques. 
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I.  Introduction 


The  Global  Positioning  System  (GPS)  is  a  satellite-based,  worldwide,  all-weather  nav¬ 
igation  and  timing  system  [1],  The  ever-increasing  reliance  on  GPS  for  navigation  and 
guidance  has  created  a  growing  awareness  of  the  need  for  adequate  protection  against 
both  unintentional  and  intentional  interference.  Jamming  is  a  procedure  that  attempts 
to  block  reception  of  the  desired  signal  by  the  intended  receiver.  In  general  terms,  it  is 
high  power  signal  that  occupies  the  same  frequency  as  the  desired  signal,  making  reception 
by  the  intended  receiver  difficult  or  impossible.  Designers  of  military  as  well  as  commer¬ 
cial  communication  systems  have,  through  the  years,  developed  numerous  anti-jamming 
techniques  to  counter  these  threats.  As  these  techniques  become  effective  for  interference 
removal  and  mitigation,  jammers  themselves  have  become  smarter  and  more  sophisticated, 
and  generate  signals,  which  are  difficult  to  combat. 

The  GPS  system  employs  BPSK-modulated  direct  sequence  spread  spectrum  (DSSS) 
signals.  The  DSSS  systems  are  implicitly  able  to  provide  a  certain  degree  of  protection 
against  intentional  or  non-intentional  jammers.  However,  in  many  cases,  the  jammer  may 
be  much  stronger  than  the  GPS  signal,  and  the  spreading  gain  might  be  insufficient  to 
decode  the  useful  data  reliably  [2].  There  are  several  methods  that  have  been  proposed  for 
interference  suppression  in  DSSS  communications  [3],  [4],  [5].  The  recent  development  of 
the  bilinear  time-frequency  distributions  (TFDs)  for  improved  signal  power  localization  in 
the  time-frequency  plane  has  motivated  several  new  effective  approaches,  based  on  instan¬ 
taneous  frequency  (IF)  estimation,  for  non-stationary  interference  excisions  [6].  One  of 
the  important  IF-based  interference  rejection  techniques  uses  the  jammer  IF  to  construct 
a  time- varying  excision  notch  filter  that  effectively  removes  the  interference  [7].  However, 
this  notch  filtering  excision  technique  causes  significant  distortions  to  the  desired  signal, 
leading  to  undesired  receiver  performance. 
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Recently,  subspacc'  projection  techniques,  which  are  also  based  on  IF  estimation,  have 
been  devised  for  non-stationary  FM  interference  excision  in  DSSS  communications  [8j.  The 
techniques  assume  clear  Jammer  time-frequency  signatures  and  rely  on  the  distinct  differ¬ 
ences  in  the  localization  properties  between  the  jammer  and  the  spread  spectrum  signals. 
The  jammer  instantaneous  frequency,  whether  provided  by  the  time-frequency  distribu¬ 
tions  or  any  other  IF  estimator,  is  used  to  form  an  interference  subspace.  Projection  can 
then  be  performed  to  excise  the  jammer  from  the  incoming  signal  prior  to  correlation  with 
the  receiver  PN  sequence.  The  result  is  improved  receiver  SINK  and  reduced  BERs. 

In  this  report,  we  apply  the  subspace  projection  techniques  as  a  pre-correlation  signal 
processing  method  to  the  FM  interference  suppression  in  GPS  receivers.  The  GPS  receiver 
and  signal  structure  impose  new  constraints  on  the  problem  since  the  spreading  code  from 
each  satellite  is  known  and  periodic  within  one  navigation  data  symbol.  This  structure  and 
the  signal  model  are  reviewed  in  Section  2.  In  Section  3,  we  depict  the  received  GPS  signal 
properties  in  time-frequency  domain.  The  SINR  of  the  GPS  receiver  implementing  the 
subspace  projection  techniques  is  derived  in  Section  4,  which  shows  improved  performance 
in  strong  interference  environments. 

II.  Signal  Model 

GPS  employs  BPSK-modulated  DSSS  signals.  The  navigation  data  is  transmitted  at  a 
symbol  rate  of  50  bps.  It  is  spread  by  a  coarse  acquisition  (C/A)  code  and  a  precision  (P) 
code.  The  C/A  code  is  a  Gold  sequence  with  a  chip  rate  of  1.023  MHz  and  a  period  of  1023 
chips,  i.e.  its  period  is  1ms,  and  there  are  20  periods  within  one  data  symbol.  The  P  code 
is  a  pseudorandom  code  at  the  rate  of  10.23  MHz  and  with  a  period  of  1  week.  These  two 
spreading  codes  are  multiplexed  in  quadrature  phases.  Figure  1  shows  the  signal  structure. 
The  carrier  LI  is  modulated  by  both  C/A  code  and  P  code,  whereas  the  carrier  L2  is  only 
modulated  by  P  code.  In  this  report,  we  will  mainly  address  the  problem  of  anti-jamming 
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for  the  C/A  code,  for  which  the  peak  power  spectral  density  exceeds  that  of  the  P  code 
by  about  13  dB  [1].  The  transmitted  GPS  signal  is  also  very  weak  with  Jammer- to-Signal 
Ratio  (JSR)  often  larger  than  40  dB  and  Signal-to-Noise  Ratio  (SNR)  in  the  range  -14  to 
-20  dB  [2],  [9].  Due  to  the  high  JSR,  the  FM  jammer  often  has  a  clear  signature  in  the 
time-frequency  domain  as  shown  in  Section  3.  As  the  P  code  is  very  weak  compared  to 
the  C/A  code,  noise  and  jammer,  we  can  ignore  its  presence  in  our  analysis. 


Fig.  1.  The  GPS  signal  structure. 


The  BPSK-modulated  DSSS  signal  may  be  expressedias  - 

=  >  7iG{-l,l}Vz  ;  Ml) 

lere  /j  represents  the  binary  information  sequence  and  Tj,  is  the  bit  interval,  which  is 
20ms  in  the  case  of  GPS  system.  The  binary  information  bit,  bi{t)  is  further  decomposed 
as  a  superposition  of  L  spreading  codes,  p(n),  pulse  shaped  by  a  unit-energy  function,  q{t), 
of  duration  of  Tc,  which  is  1/1023  ms  in  the  case  of  C/A  code.  Accordingly, 

^i(0  =  S  P(")  ~  ’^^c)  (2) 

n=l 

The  signal  for  one  data  bit  at  the  receiver,  after  demodulation,  and  sampling  at  chip 
rate,  becomes 

x(n)  =  p(n)  -F  u;(n)  +j(n)  1  <  n  <  L  (3) 
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where  piv)  is  the  chip  sequence.  u'{7i)  is  the  white  noise,  and  j{n)  is  the  interfeiing  signal. 
The  above  equation  can  b('  written  in  the  vector  form 

X  =  p  +  w  +  j  (4) 


w'here 


X 

P 

w 

j 


x(l)  .t(2)  x(31  x{L) 

p{l)  p{2)  p(3)  p[L) 


T 

T 


U-(1)  w{2)  w(3)  v'iL) 

I'l 


3(1)  ./(2)  3(3)  j{L) 


All  vectors  are  of  dimension  L  x  1,  and  ‘T  denotes  vector  or  matrix  transposition.  It 
should  be  noted  that  the  P  vector  is  real,  whereas  all  other  vectors  in  the  above  equation 


have  complex  enterics. 


III.  Periodic  Signal  Plus  Jammer  in  the  Time-Frequency  Domain 

For  GPS  C/A  code,  the  PN  sequence  is  periodic.  The  PN  code  of  length  1023  repeats 
itself  20  times  within  one  symbol  of  the  50  bps  navigation  data.  Consequently  it  is  no 
longer  of  a  continuous  spectrum  in  the  frequency  domain,  but  rather  of  spectral  lines. 
The  case  is  the  same  for  periodic  jammers.  Figure  2  and  Figure  3  show  the  effect  of 
periodicit}-  of  the  signal  and  the  jammer  on  their  respective  power  distribution  over  time 
and  frequency  using  Wigner-^'ille  distribution.  In  both  figures,  a  PN  sequence  of  length 
32  samples  that  repeats  8  times  is  used.  A  non-periodic  chirp  jammer  of  a  50dB  JSR 
(jammer-to-signal  ratio)  is  added  in  Figure  2.  A  periodic  chirp  jammer  of  50  dB  JSR 
with  the  same  period  as  the  C/A  code  is  included  in  Figure  3.  We  note  that  the  chosen 
value  of  50dB  JSR  has  a  practical  significance.  The  spread  spectrum  systems  in  a  typical 
GPS  C/A  code  receiver  can  tolerate  a  narrowband  interference  of  approximately  40  dB 
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JSR  without  interference  mitigation  processing.  However,  field  tests  show  that  jammer 
strength  often  exceeds  that  number  due  to  the  weakness  of  the  signal.  SNR  in  both  figures 
are  -20dB.  which  is  also  close  to  its  practical  value  [2],  [9].  Due  to  high  JSR,  the  jammer 
is  dominant  in  both  figures.  From  Figure  3,  it  is  clear  that  the  periodicity  of  the  jammer 
brings  more  difficulty  to  IF  estimation  than  the  non-periodic  jammers.  This  problem  can 
be  solved  by  applying  a  short  data  window  when  using  Wigner-Ville  distribution.  Note 
that  the  window  length  should  be  less  than  the  jammer  period.  Figure  4  shows  the  result 
of  applying  a  window  of  length  31  to  the  same  data  used  in  Fig.  3.  It  is  evident  from  the 
Fig.  4  that  the  horizontal  discrete  harmonic  lines  have  disappeared. 


contour  of  periodic  PN  +  nonperiodic  chirp  jammer 


Fig.  2.  Periodic  signal  corrupted  by  a  non-periodic  jaininer  in  time-frequency  domain 


IV.  GPS  Anti-Jamming  Using  Projection  Techniques 

The  concept  of  subspace  projection  for  instantaneously  narrowband  jammer  suppression 
is  to  remove  the  jammer  components  from  the  received  data  by  projecting  it  onto  the 
subspace  that  is  orthogonal  to  the  jammer  subspace,  as  illustrated  in  Fig.  5. 

Once  the  instantaneous  frequency  (IF)  of  the  non-stationary  jammer  is  estimated  from 
the  time-frequency  domain,  or  by  using  any  other  IF  estimator  [10],  [11],  [12],  [13],  the 
interference  signal  vector  j  in  (4)  can  be  constructed,  up  to  ambiguity  in  phase  and  possibly 


9 


contour  of  periodic  PN  +  periodic  chirp  jammer 


“f 


60  100  160  200  260 
time  domain 


Fig.  3.  Periodic  signal  corrupted  by  a  periodic  jaininer  in  time-frequency  domain 


cor^touf  of  periodic  PN  4  periodic  chirp  jammer,  with  window 


Fig.  4.  Periodic  signal  corrupted  by  a  periodic  jainnier  in  time-frequency  domain  (with  window) 


in  amplitude.  In  the  proposed  interference  excision  approach,  the  data  vector  is  partitioned 
into  Q  blocks,  each  of  length  P,  i.e.  L=PQ.  For  the  GPS  C/A  code,  Q=20,  P=1023,  and  all 
Q  blocks  are  identical,  i.e.,  the  signal  PN  sequence  is  periodic.  Block-processing  provides 
the  flexibility  to  discard  the  portions  of  the  data  bit,  over  which  there  are  significant 
errors  in  the  IF  estimates.  The  orthogonal  projection  method  makes  use  of  the  fact  that, 
in  each  block,  the  jammer  has  a  one-dimensional  subspace  J  in  the  P-dimensional  space 
V,  which  is  spanned  by  the  received  data  vector.  The  interference  can  be  removed  from 
each  block  by  projecting  the  received  data  on  the  corresponding  orthogonal  subspace  Q  of 
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Signal  Subspace 


i  Orthogonal  Subspace  G 

V  Px 


/ 

Fig.  5.  Jairiiiier  excision  by  subspace  projection 


interference  subspace  J.  The  subspace  J  is  estimated  using  the  IF  information.  The 
projection  matrix  for  the  block  is  given  by 

Vfc  =  I  -  Ufc  uf  (5) 


The  vector  Ujt  is  the  unit  norm  basis  vector  in  the  direction  of  the  interference  vector  of  the 
block,  and  ‘H’  denotes  vector  or  matrix  Hermitian.  Since  the  FM  jammer  signals  are 
uniquely  characterized  by  their  IFs,  the  f**  FM  jammer  in  the  block  can  be  expressed 
as 


Ukii)  =  -^^expljcjfkii)] 

The  result  of  the  projection  over  the  k^^  data  block  is 


(6) 


Xfc  =  V*:  Xfc 


(7) 


where  x*,  is  the  input  data  vector.  Using  the  three  different  components  that  make  up  the 
input  vector  in  (4),  the  output  of  the  projection  filter  V*,  can  be  written  as 

Xfc  =  Vfc  [pfc  +  Wfc+jJ  (8) 

The  noise  is  assumed  to  be  complex  white  Gaussian  with  zero-mean, 

E[w{n)]  —  0,  E[w{n)*w{n  + 1)]  =  (9) 


n 


Since  we  assume  total  interference  excision  through  the  projection  operation,  then 

=  0,  Xa- =  Vfc  p*.  +  V^.  w*.  (10) 

The  decision  variable  is  the  real  part  of  y  that  is  obtained  by  correlating  the  filter 

output  xa:  with  the  corresponding  block  of  the  receiver  PN  sequence  and  summing  the 
results  over  the  K  blocks.  That  is. 

y  =  (11) 

k-O 

Since  the  PN  code  is  periodic,  we  can  strip  off  the  subscrijjt  k  in  pk-  The  above  variable 
can  be  written  in  terms  of  the  constituent  signals  as 

y  =  X!  P^Vfcp  +  51  w^Va-p  a  yi  +  y-,  (12) 

k=0  k=0 

where  yi  and  y2  are  the  contributions  of  the  PN  and  noise  sequences  to  the  decision 

variable,  respectively.  In  [8],  yi  is  considered  as  a  random  variable.  However,  in  GPS 

system,  due  to  the  fact  that  each  satellite  is  assigned  a  fixed  Gold  code  [1],  and  that  the 
Gold  code  is  the  same  for  every  navigation  data  symbol,  yi  can  no  longer  be  treated  as 
a  random  variable,  but  rather  a  deterministic  value.  This  is  a  key  difference  between  the 
GPS  system  and  other  spread  spectrum  systems.  The  value  of  yi  is  given  by 


Define 


Q-i 

yi  ^  Y.  P^VaP 

A:=0 

=  5Z  P^(I  -  Ufcuf  )p 

k=0 

=  5I(p^p  “ 

*:=0 

Q-1 

=  QP  -  YiP^^k-iikP) 

k=0 


3k  = 


P^ua 

^/P 


(13) 


(14) 
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as  the  correlation  coefficient  between  the  PN  sequence  vector  p  and  the  jammer  vector  u. 
Pk  reflects  the  the  component  of  the  signal  that  is  in  the  jammer  subspace,  and  represents 
the  degree  of  resemblance  between  the  signal  sequence  and  the  jammer  sequence.  Since 
the  signal  is  a  PN  sequence,  and  the  jammer  is  a  non-stationary  FM  signal,  the  correlation 
coefficient  is  typically  very  small.  With  the  above  definition,  j/i  can  be  expressed  as 

EIAP)  (15) 

*:=0 

From  (15),  it  is  clear  that  j/i  is  a  real  value,  which  is  the  result  of  the  fact  that  the  projection 
matrix  is  Hermitian.  With  the  assumptions  in  (9).  y2  is  complex  white  Gaussian  with 
zero-mean.  Therefore, 


4  =  -^[1^21"] 

=  E  (5Z  w^Vfcp)"(^  w^V,p) 

*=0  1=0 

=  YlYl  P^^kE  [wfcwf  1  V(P 
k=0  1=0 

=  Yi  P^^kE  [wfcwf  ]  VkP 

k=0 

A:=0 

=  a^Y  P^^kP  ^  CT^yi  (16) 

fc=0 

the  above  equations  make  use  of  the  noise  assumptions  in  (9)  and  the  properties  of  the 
projection  matrix.  The  decision  variable  is  the  real  part  of  y.  Consequently,  y^  is  given 
by 

yr  =  yi  -t-  Re{y2}  (17) 


where  Re{y2}  denotes  the  real  part  of  y2.  i?e{y2}  is  real  white  Gaussian  with  zero-mean 


and  variance  .  Therefore,  the  SINK  is 


SINR  = 


var{Re{y2}} 
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(18) 


2P{Q  -  Etd  IAI-) 

cr'^ 

In  the  absence  of  jammers,  no  excision  is  necessary,  and  the  SINR(SNR)  of  the  receiver 
output  will  become  2PQjo^,  which  represents  the  upper  bound  for  the  anti-jamming  per¬ 
formance.  Clearly,  —  is  the  reduction  in  the  receiver  performance  caused  by  the 

proposed  jammer  suppression  techniques.  It  reflects  the  energy  of  the  power  of  the  signal 
component  that  is  in  the  jammer  subspace.  If  the  jammer  and  spread  spectrum  signals 
are  orthogonal,  i.e.,  their  correlation  coefficient  l/l]  ==  0.  then  interference  suppression  is 
achieved  with  no  loss  in  performance.  However,  as  stated  above,  in  the  general  case,  is 
often  very  small,  so  the  projection  technique  can  excise  FM  jammers  effectively  with  only 
very  insignificant  signal  loss.  The  lower  bound  of  SINR  is  zero  and  corresponds  to  |/3|  =  1. 
This  case  requires  the  jammer  to  assume  the  C/A  code,  i.e.,  identical  and  synchronous 
with  actual  one.  Figure  6  depicts  the  theoretical  SINR  in  (18),  its  upper  bound,  and 
estimated  values  using  computer  simulation.  The  SNR  assumes  five  different  values  [-25, 
-20,  -15,  -10,-5]  dB.  In  this  figure,  the  signal  is  the  Gold  code  of  satellite  SV#1.  and  the 
jammer  is  a  periodic  chirp  FM  signal  with  frequency  0-0.5  and  has  the  same  period  as 
the  C/A  code.  For  this  case,  the  correlation  coefficient  /?  is  very  small,  |/1|  =  0.0387.  JSR 
used  in  the  computer  simulation  is  set  to  50dB.  Due  to  the  large  computation  involved, 
we  have  used  1000  realizations  for  each  SNR  value.  Figure  6  demonstrates  that  the  theo¬ 
retical  value  of  SINRs  is  almost  the  same  as  the  upper  bound  and  both  are  very  close  to 
the  simulation  result.  In  the  simulation  as  well  as  in  the  derivation  of  equation  (18),  we 
have  assumed  exact  knowledge  of  the  jammer  IF.  Inaccuracies  in  the  IF  estimation  will 
have  an  effect  on  the  receiver  performance  [8]. 
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Fig.  6.  Receiver  SINK  vs  SNR. 


V.  Conclusions 

GPS  receivers  are  vulnerable  to  strong  interferences.  In  this  report,  subspace  projection 
techniques  are  adapted  for  the  anti-FM  jamming  GPS  receiver.  These  techniques  are  based 
on  IF  estimation  of  the  jammer  signal,  which  can  be  easily  achieved,  providing  that  the 
C/A  code  and  the  jammer  have  distinct  time-frequency  signatures.  The  IF  information  is 
used  to  construct  the  FM  interference  subspace  which,  because  of  signal  nonstationarities, 
is  otherwise  difficult  to  obtain.  Due  to  the  characteristic  of  the  GPS  spread  spectrum  signal 
structure  and  the  fact  that  the  C/A  codes  are  fixed  for  the  different  satellites  and  known 
to  all,  the  analysis  of  the  receiver  SINK  becomes  different  from  common  spread  spectrum 
systems.  The  theoretical  and  simulation  results  suggest  that  the  subspace  projection 
techniques  can  effectively  excise  FM  jammers  for  GPS  receivers  with  insignificant  loss  in 
the  spreading  gain. 
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Chapter  2 


Mitigation  of  Periodic  Interferers  in  GPS  Receivers  Using 
Subspace  Projection  Techniques 


Abstract 

Frequency  modulated  signals  in  the  frequency  band  1.217-1.238  GHz  and  1.565- 
1.586  GHz  present  a  source  of  interference  to  the  GPS,  which  should  be  properly 
mitigated.  In  this  report,  we  derive  the  Signal-to-Interference-and-Noise  Ratio 
(SINK)  of  the  GPS  receiver  implementing  subspace  projection  techniques  for 
suppression  of  FM  jammers.  We  consider  the  general  case  in  which  the  jammer  may 
have  equal  or  different  cycles  than  the  coarse  acquisition  (C/A)  code  of  the  GPS 
signals.  It  is  shown  that  the  weak  correlations  between  the  FM  interference  and  the 
Gold  codes  allow  effective  interference  cancellation  without  significant  loss  of  the 
desired  signal 
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I.  Introduction 


Subspace  projection  techniques  based  on  time-frequency  distributions  have  been 
employed  for  suppression  of  non-stationary  FM  interference  in  broadband 
communication  platforms  [l]-[4].  Most  recently,  they  have  been  applied  to  Global 
Positioning  System  (GPS)  with  single  [5]  and  multi-sensor  [6]  receivers.  These 
techniques  assume  clear  jammer  time-frequency  signatures  and  rely  on  the  distinct 
differences  in  the  localization  properties  between  the  interference  FM  waveforms  and  the 
coarse  acquisition  (C/A)  Gold  codes  of  the  GPS  signals.  The  FM  jammer  instantaneous 
frequency  (IF),  whether  provided  by  the  time-frequency  distributions  or  any  other  IF 
estimator  [7]  [8],  is  used  to  define  the  temporal  signature  of  the  interference,  which  is  in 
turn  used  to  construct  the  interference  subspace.  The  respective  projection  matrix  is  used 
to  excise  the  jammer  power  in  the  incoming  signal  prior  to  correlation  with  the  receiver 
C/A  codes.  The  result  is  improved  receiver  signal-to-interference-plus-noise  ratio 
(SINK)  and  reduced  BERs. 

In  this  report,  we  generalize  the  results  in  Chapter  1  and  reference  [5]  by  considering 
jammer  signals  with  different  periodic  structures  from  that  of  the  GPS  C/A  codes.  In  the 
underlying  problem,  we  deal  with  the  case  in  which  multiplicities  of  the  jammer  period 
span  a  finite  number  of  the  GPS  information  symbols.  Therefore,  unlike  the  previous 
work  that  assumes  periodic  synchronization  between  the  interference  and  the  desired 
signal,  the  generalization  herein  allows  different  portions  of  the  jammer  signal  to  infringe 
on  different  symbols  of  the  GPS  C/A  code.  It  is  shown,  however,  that  due  to  the  weak 
correlation  between  the  FM  waveforms  and  the  Gold  codes,  the  GPS  receiver 
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implementing  subspace  projections  is  robust  to  FM  jammer  periodic  patterns  and 
repetition  cycles,  achieving  full  interference  suppression  with  no  significant  performance 
degradation  from  that  obtained  in  the  optimum  interference  free  environment. 

II.  SINK  Analysis  for  Periodic  Jammers 

The  derivation  in  Chapter  1  implicitly  assumes  that  the  jammer  period,  7j,  is  equal  to  the 
symbol  length,  T^,,  of  the  GPS  signal,  i.e.,  Tj=Tf,=PQ.  This  assumption  limits  the  SINK 
result 

(  Q-\  W 

SINR  =  2P  (T-  (1) 

I  A=o  j/ 

to  the  special  case  where  the  jammer  and  GPS  signals  have  the  same  cycle.  Considering 
the  more  general  case,  the  jammer  is  presumed  to  be  a  periodic  signal  with  Tj=-^PQ, 

where  N  and  M  are  integers  and  relatively  prime.  Recall  that  P=1023  is  the  period  of  the 
GPS  spreading  code,  and  Q=20  is  the  number  of  times  the  code  repeats  itself  over  one 
symbol.  From  the  above  definition,  M  jammer  periods  extend  over  N  symbols  of  the 
GPS  signal,  and,  as  such,  different  segments  of  the  jammer  signature  will  infringe  on 
different  GPS  symbols.  Therefore,  in  the  analysis  herein,  we  add  the  subscript  /  to 
associate  the  receiver  variables  with  the  k'''  symbol.  It  is  straightforward  to  show  that  the 
correlation  output  at  the  receiver  yields 

>’,■  =  SP^V,.tP  +Xw^(n)V,.^.pAy.,  +  y,.^  (2) 

i=0  k=0 

The  decision  variable  is  the  real  part  of  y,.  It  can  be  shown  that 

y,  =  m (3) 

Jk=0 
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On  the  other  hand,  the  correlation  output  due  to  the  noise,  y,:,  is  a  complex  Gaussian 
zero-mean  random  variable,  and  its  variance  can  be  readily  obtained  as 


=o-'/’(!2-X|4r)  (4) 

k=0 

It  is  noted  that  since  there  are  N  symbols  for  every  M  jammer  periods,  and  M^N,  both 
variables  y,)  and  y,2  assume  different  values  over  N  consecutive  GPS  symbols.  The 
jammer  can  then  be  cast  as  symbol-dependent,  assuming  N  distinct  waveforms.  In  this 
case,  one  simple  measure  of  the  receiver  performance  is  to  average  the  SINR  over  N 


consecutive  symbols,  i.e.. 


SINR^^,  =  E[SINR,] 


=  ^P,(SINR,\J,)P^(J.)  (5) 

/=i 

1  'V 

=— y  s/m 

where  SINRi  and  SINRav  denote,  respectively,  the  receiver  signal-to-interference-and- 
noise  ratio  over  the  i''‘  symbol  and  the  average  receiver  SINR.  In  the  above  equation, 
SINRj  is  treated  as  a  discrete  random  variable  that  takes  N  possible  values  with  equal 

probability.  7,(;  =  i . A^)are  the  segments  of  the  jammer  signal  over  N  consecutive 

symbols.  In  equation  (5)  Pr(x)  denotes  the  probability  of  the  event  x,  and  Pr(y/j=l/A^.  The 


SINRi  is 


2/’(e-2lA.r) 


SINR.  - - ^ 

<7 


Accordingly, 


S/m  „  = 


iV  ,=1 
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The  above  expression,  although  simple  to  calculate,  smoothes  out  high  and  low  SINR 
values.  In  this  regard,  the  average  value  in  (7)  does  not  properly  penalize  poor  or  reward 
good  receiver  performance.  Further,  it  is  difficult  to  establish  a  relationship  between  the 
receiver  SINRav  and  its  BER.  Most  importantly,  expression  (7)  does  not  account  for  the 
self-noise  term  that  reflects  the  level  of  signal  distortion  produced  by  the  induced 
correlation  of  the  code  chips  as  a  result  of  the  excision  process.  Hence,  a  more  proper 
way  to  measure  the  receiver  performance  is  to  deal  with  v/  as  a  random  variable.  In  this 
case  the  average  receiver  signal  to  interference  plus  noise  ratio  is  referred  to  as  siNR  to 
distinguish  it  from  equation  (1).  We  assume  that  symbol  “1  is  transmitted  and 
contaminated  by  one  of  N  possible  jammer  signals  occurring  with  the  same  probability. 
In  this  case,  the  mean  value  and  the  variance  of  the  correlator  output  due  to  the  GPS 
signal  can  be  derived  as 


E[y\]  =  I  J,]PAJ,)  =  I 

iV  1=1 


/V 


*=0 


N 


1  V 


<T-  =£[)fl-£=(.v,J 

4ip=(e-SIAr)=-;>'e=(i-^il;iAr)= 

ly  ,=i  *=o  1=1  *=o 

d3  \  C-l  p2  N  e-l 

=^S(Si«.r)=-4<ssiAr)^ 

ly  /=i  *=o  1=1  *=o 

Similarly,  the  average  noise  power  is 


(8) 


(9) 
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<  =E[yl]  =  ^E[y;\J^]Pr{Ji) 


/=1 


Q-l 


2;o-=p(e-ElA.r) 


i=l 


it=0 


N 


From  the  above  and  ( 1 )  the  average  SINR  is  given  by. 


(10) 


SINR  = 


E\y,] 

a^.  +-o^. 

■M  2  •'* 


1  w  2-1 

1V^,=:U=0 


This  expression  represents  the  SINR  of  the  receiver  implementing  subspace 
projections  to  remove  a  periodic  jammer  and  is  also  valid  for  the  case  in  which  the 
jammer  assumes  N  possible  waveforms.  In  the  case  that  the  jammer  has  the  same  period 
as  the  GPS  data  symbol,  then  N=l  and  the  above  equation  simplifies  to  single  antenna 
case.  Comparing  (11)  to  (7)  and  SINR  expression  in  Chapter  1,  it  is  clear  that  (11) 
includes  the  self-noise  component  o-*  that  arises  due  to  the  differences  in  the  distortion 
effects  of  interference  excision  on  the  GPS  signal  over  the  N  symbols.  In  the  absence  of 
jamming,  no  excision  is  necessary,  and  the  SINR  of  the  receiver  output  becomes  2PQ/<y, 
which  represents  the  upper  bound  of  the  interference  suppression  performance. 
Moreover,  if  the  jammer  and  the  spreading  codes  are  orthogonal,  i.e.,  the 

interference  suppression  is  also  achieved  with  no  loss  in  optimum  receiver  performance. 


23 


It  is  noted,  however,  that  the  values  of  the  cross-correlation  coefficient.  |/?,*| .  between  the 
PN  sequence  signal  and  the  non-stalionary  FM  jammer  are  typically  very  small.  This 
allows  the  proposed  projection  technique  to  excise  FM  jammers  effectively  with 
insignificant  signal  loss.  Computer  simulations  show  that  |/?,^|  ranges  from  0  to  0.14. 

With  the.se  values,  the  self-noise  <77^  is  negligible  compared  to  the  Gaussian  noise  for  the 
low  SNR  conditions  that  often  prevail  in  GPS  environment.  In  this  case,  equation  (11) 
can  be  simplified  to  the  following 


SINR 


g-[yi] 

0 


2PQ 


1  N  0-1  .  ^ 


NQ..U-0 


(12) 


which  is  similar  to  the  SINR  expression  in  Chapter  1  and  has  the  same  form  as  (7). 


Therefore,  SIMR^v  and  SINR  approximately  yield  the  same  performance  measure. 


III.  Simulation  Results 

Fig.  1  plots  the  receiver  SINR  vs  SNR  according  to  (11)  for  the  two  cases  of  N=l  and 
N/M=5/3.  In  both  cases,  the  normalized  start  and  end  frequencies  of  the  chirp  jammer  are 
0  and  0.5,  respectively.  The  SNR  values  range  from  -25dB  to  -5dB,  and  the  GPS  signal 
is  the  Gold  code  of  satellite  SV#1.  It  is  clear  that  the  period  of  the  jammer  has  little 
effects  on  the  result  of  interference  suppression  performance,  as  both  SINR  curves  are 
very  close  to  the  upper  bound.  From  Fig.  1 ,  we  can  also  observe  that  the  SINR  change 
linearly  with  the  input  SNR,  which  can  be  easily  recognized  from  (12).  Fig.  2  shows  the 
values  for  the  underlying  example.  It  is  evident  from  this  figure  that  there  is  no  clear 

variation  patterns  of  the  cross-correlation  coefficients.  The  range  values  of  I/?,*!  do  not 
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change  greatly  over  different  symbols.  Experiments  with  different  N,  M,  chirp  rate,  and 
satellite  signals  have  given  similar  results. 

IV.  Conclusions 

Subspace  projection  techniques  were  employed  as  a  pre-correlation  signal  processing 
method  in  the  GPS  receiver  to  excise  non-stationary  FM  jammers.  Previous 
developments  in  receiver  performance  evaluation  for  periodic  interferers  have  only 
considered  the  special  case  of  jammers  having  the  same  period  as  the  GPS  symbol.  This 
report  establishes  the  derivation  for  the  receiver  SINR  expression  for  non-stationary  FM 
jammers  with  general  periodic  structures.  It  was  argued  that  in  the  general  case  the  FM 
jammer  contaminates  each  GPS  symbol  with  a  finite  number  of  possible  waveforms. 
Two  averaging  measures  to  evaluate  the  overall  SINK  performance  have  been  discussed 
and  compared.  The  simulations  included  in  this  report  employed  chirp  interferers  and  it 
was  shown  that  these  two  measures  give  similar  performance  for  this  jammer  class.  The 
results  suggest  that  the  subspace  projection  method  can  excise  periodic  FM  jammers 
effectively  with  insignificant  losses  to  the  desired  signal. 
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Receiver  output  SIN  R  after  subsjgace  projection  (dB) 


SNR(dB) 


Figure  1.  Output  SINK  vs  SNR. 
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Chapter  3 


Subspace  Array  Processing  for  the  Suppression  of 
FM  Jammers  in  GPS  Receivers 


Abstract 

This  report  applies  the  subspace  projection  array  processing  techniques  for  suppression 
of  frequency  modulated  (FM)  jammers  in  GPS  receivers.  The  FM  jammers  are 
instantaneous  narrowband  and  have  clear  time-frequency  (t-f)  signatures  that  are  distinct 
from  the  GPS  C/A  spread  spectrum  code.  In  the  proposed  technique,  the  instantaneous 
frequency  (IF)  of  the  jammer  is  estimated  and  used  to  construct  the  jammer  subspace. 
With  a  multi-sensor  receiver,  both  spatial  and  time-frequency  signatures  of  signal  arrivals 
are  used  for  effective  interference  suppression.  This  report  considers  the  deterministic 
nature  of  the  GPS  C/A  code.  We  derive  the  receiver  SINR,  which  shows  improved 
performance  in  strong  interference  environments. 
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I.  Introduction 


Compared  with  the  subspace  projection  techniques  in  the  single-sensor  case,  the  use  of 
multi-sensor  array  greatly  increases  the  dimension  of  the  available  signal  subspace.  It 
allows  both  the  distinctions  in  the  spatial  and  temporal  signatures  of  the  GPS  signals  from 
those  of  the  interferers  to  play  equal  roles  in  suppressing  the  jammer  with  a  minimum 
distortion  of  the  desired  signal.  In  this  report,  we  examine  the  applicability  of  multi¬ 
sensor  subspace  projection  techniques  for  suppressing  nonstationary  jammers  in  GPS 
receivers. 

We  mainly  address  the  problem  of  anti-jamming  for  the  C/A  code.  This  code  has 
interesting  properties  due  to  the  fact  that  it  is  periodic  within  one  data  symbol  and  fixed 
for  each  satellite  signal.  It  is  shown  that  this  periodicity  has  a  considerable  effect  on  the 
receiver  performance. 

II.  Subspace  projection  array  processing 

In  GPS,  the  PN  sequence  of  length  L  (1023)  repeats  itself  Q  (20)  times  within  one 
symbol  of  the  50  bps  navigation  data.  We  use  discrete-time  form,  where  all  the  signals 
are  sampled  at  the  chip-rate  of  the  C/A  code.  We  consider  an  antenna  array  of  N  sensors, 
and  the  communication  channel  is  restricted  to  flat-fading.  In  the  proposed  interference 
excision  approach,  the  LNQ  sensor  output  samples  are  partitioned  into  Q  blocks,  each  of 
L  chips  and  LN  samples.  The  jammer  can  be  consecutively  removed  from  the  20  blocks 
that  constitute  on  symbol.  This  is  achieved  by  projecting  the  received  data  in  each  block 
on  the  corresponding  orthogonal  subspace  of  the  jammer.  The  jammer-free  signal  is  then 
correlated  with  the  replica  PN  sequence  on  a  symbol-by-symbol  basis.  We  first  consider 
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the  subspace  projection  within  each  block.  The  array  output  vector  at  the  k-th  sample  is 
given  by 


x(k)  =  x^(k)  +  x^^ik)  +  b{k) 

^  (1) 
=  c{k)h  +  ^AiUjik)ai  +  b{k) 
i=\ 

where  x,,  and  b  are  the  signal,  the  jammer  and  the  white  Gaussian  noise 
contributions,  respectively,  h  is  the  signal  spatial  signature,  and  c(k)  is  the  spreading  PN 
sequence.  The  number  of  jammers  is  U.  All  jammers  are  considered  as  instantaneously 
narrowband  FM  signals  with  constant  amplitude  u^{k)=txp[j(p.{k)].  A.  and  a,  are  the  i-th 
jammer  amplitude  and  the  spatial  signature,  respectively.  Furthermore,  we  normalize  the 
channels  and  set  ||h||^  =  N  and  ||a|l^  =  N  ,  where  ||^is  the  Frobenius  norm  of  a  vector.  The 
noise  vector  b(A:)  is  zero-mean,  temporally  and  spatially  white  with 

E[b{k)b'^{k  +  l)]  =  0  ^2) 

E[E[b{k)b^  (k  +  /)]  =  c7^^(/)/yv 

where  is  the  noise  power,  and  In  is  the  NxN  identity  matrix.  Using  L  sequential 
array  vector  samples  within  the  block,  we  obtain  the  following  LA-by-1  vector 

X=[x^(l)  x^(2)  ...  x^iL)f 

=  X^  +  X,^+B  ^2) 

The  vector  X„  consists  of  the  U  jammer  signals,  and  is  expressed  as 


x„=I,A.yi 


Vi=[u,(l)  «,(2)  ...  u^m'  ®a, 
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where  ®  denotes  the  Kronecker  product.  Therefore, 

y=[V,  V2-^t/]  '  (6) 

spans  the  jammer  signal  subspace,  and  its  orthogonal  subspace  projection  matrix  is 
given  by 

p=hN-y(y^yy^y''  =hN—^yy^  o) 

LN 

The  projection  of  the  received  signal  vector  onto  the  orthogonal  subspace  yields 

X^=PX=  PX,  +  PB  (8) 

which  excises  the  jammers  entirely.  The  signal  vector  can  be  rewritten  as 

X,  =[c(l)  c(2)  ...  c(L)f  ®h  A  q  (9) 

where  the  vector  q  represents  the  spatial-temporal  signature  of  the  GPS  signal.  The 
result  of  despreading  in  the  subspace  projection  based  array  system  over  one  block  is 

y  =  q^X^=q^Pq  +  q^PB  A  +  y2  (1®) 

where  yi  and  y2  are  the  contributions  of  the  PN  and  the  noise  sequences  to  the  decision 
variable,  respectively.  For  simplification,  we  assume  that  the  jammers  share  the  same 
period  as  the  GPS  data  symbol.  In  GPS  systems,  due  to  the  fact  that  each  satellite  is 
assigned  a  fixed  Gold  code  [3],  and  that  the  Gold  code  is  the  same  for  every  navigation 
data  symbol,  yi  is  a  deterministic  value,  rather  than  a  random  variable  as  it  is  the  case  in 
many  spread  spectrum  applications.  The  value  of  yi  is  given  by 
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where 


Define 


LN 


=  V,  ...  v„] 


v, 

V. 

yu 


=  LN-—q 
LN 


H 


fu 


Sv,v« 


/=1 


=  LN - —Yg'^ViVi^q 


=  (c®hfiUi  <8>a,)  =  (c^u;)(h^a,) 


«/=■ 


N 


(11) 


(12) 


(13) 


as  the  spatial  cross-correlation  coefficient  between  the  signal  and  the  i-th  jammer,  and 


T 

L 


(14) 


as  the  temporal  cross-correlation  coefficient  between  the  PN  sequence  and  the  i-th 
jammer  vector.  Therefore, 

q^Vi=LNa,P,  (15) 

JL.  ... 

(16) 

From  the  above  equation,  it  is  clear  that  yi  is  real,  which  is  due  to  the  Hermitian 
property  of  the  projection  matrix  P.  With  the  assumptions  in  (2),  y2  is  complex  Gaussian 
with  zero-mean.  Therefore, 


u 

.v,=L7v-LA^XKrlAr 

^LN 

1-XKI  IAI 

(=1 

1  1=1  j 
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£[v]=  Vi  =LN 


f  u  .  \ 

lAl 


(17) 


Var[y]  =  Var[y2]  =  E[\y2f] 

=  E[q^PBB^P^g]  =  q^PE[BB"]Pq 
=  a^q^PPq  =  CT^q^  Pq 


=  a-y^  =(J-LN 


(18) 
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1 

-SKI  1/ 
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/=i 
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The  above  two  equations  are  derived  for  only  one  block  of  the  signal  symbol.  Below, 
we  add  subscript  m  to  identify  y  with  block  m  (m=l,  2,  Q).  By  summing  all  Q  blocks, 

we  obtain  the  output  of  the  symbol-level  despreading. 


Q 

y  =  X 

/«=1 


(19) 


Since  y^  are  Gaussian  with  zero-mean,  y  is  also  a  zero-mean  Gaussian  random 
variable.  The  decision  variable  y^  is  the  real  part  of  y. 


The  expected  value  of  y^  is 


y^  =  Re[y] 


Q 


(20) 


£U.]  =  £[yl  =  2;y, 


■'m 

m=l 

f  Q  U 


=  LN 


=  LN 


2-XZK/l  \fin 
m=l  i=] 

u  (  .  Q 

Q-%  Kl 


(21) 


1=1 


m=l 


where  a^.  and  are  the  spatial  and  temporal  cross-correlation  coefficients  between 
the  signal  and  the  i-th  jammer  over  block  m.  Since  the  changes  in  the  spatial  signatures  of 
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the  signals  and  jammers  are  very  small  compared  with  the  period  of  one  block  (1  ms), 
can  be  simplified  to  cr.  The  variance  of  }\  is 

Q 

(22) 


1  t  ^ 

')  I  0  r  v"’  "• 

<7,.  =-(J~  =  -  >  (t; 

.'r  O  -  T  'm 
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t/  -  e 
|2 


Kl  XI A, 

/=!  m=\ 


Therefore,  the  receiver  SINR  expression  after  projection  and  despreading  is  given  by 


SINR^ 


E\y] 

Var{y] 


u 

f  .  Q 

2LN 

2-1 

Kl'  XlAmf 

/=1 

1  JJ 

(23) 


The  temporal  and  spatial  coefficients  appear  as  multiplicative  products  in  (23).  This 
implies  that  the  spatial  and  temporal  signatures  play  equal  roles  in  the  receiver 
performance.  In  the  absence  of  the  jammers,  no  excision  is  necessary,  and  the  SINR  of 
the  receiver  output  will  become  ILNQ/cr,  which  represents  the  upper  bound  of  the 
interference  suppression  performance.  Clearly,  the  term 


u  {  .  Q 

ZLN^  |a,| 


1=1 


m=l 


(24) 


in  equation  (23)  is  the  reduction  in  the  receiver  performance  caused  by  the  proposed 
interference  suppression  technique.  It  reflects  the  energy  of  the  signal  component  that  is 
in  the  jammer  subspace.  It  is  important  to  note  that  if  the  jammers  and  the  DSSS  signal 
are  orthogonal,  either  in  spatial  domain  (a;=0)  or  in  temporal  domain  (y5„,=0),  the 
interference  excision  is  achieved  with  no  loss  in  performance.  In  the  general  case,  are 
very  small  and  much  smaller  than  a..  Therefore,  the  difference  in  the  temporal  signatures 
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effectively  with  only  insignificant  signal  loss.  The  spatial  cross-correlation  coefficients  a. 
are  fractional  values  and,  as  such,  further  reduce  the  undesired  term  in  (24). 

Interference  suppression  using  arrays  is  improved  in  several  ways.  First,  the 
employment  of  an  antenna  array  can  lead  to  an  accurate  IF  estimation  of  the  jammers  [9]. 
Second,  in  comparison  with  the  result  of  the  single  sensor  case  [10], 


multi-sensor  receivers,  at  minimum,  improve  SINR  by  the  array  gain.  This  is  true, 
independent  of  the  underlying  fading  channels  and  scattering  environment.  As  discussed 
above,  compared  to  the  temporal  signatures,  the  contributions  of  spatial  cross-correlation 
coefficients  a,  in  influencing  the  receiver  SINR  are,  in  most  cases,  insignificant. 

III.  Simulation  Results 

In  this  section,  we  present  some  simulation  results  illustrating  the  properties  of  the 
performance  of  the  proposed  projection  technique.  Fig.  1  depicts  the  SINR  expression  as 
a  function  of  the  input  SNR.  We  consider  two  uncorrelated  chirp  jammers.  The  angle-of- 
arrival  (AOA)  of  the  satellite  signal  and  the  jammers  are  5®,  40",  and  60",  respectively.  A 
two-element  array  is  considered  with  half-wavelength  spacing.  The  satellite  PN  sequence 
is  the  Gold  code  of  satellite  SV#1,  and  the  normalized  frequency  of  the  jammers  are  from 
0.01  to  0.2  and  from  0.5  to  0.3,  respectively.  The  SINR  of  the  single  sensor  case  is  also 
plotted  for  comparison.  The  array  gain  is  evident  in  Fig.  1.  In  this  example,  Iq!;I=0.643, 
la^l=0.340,  is  in  the  range  [0.0049-0.0604],  and  is  in  the  range  [0.0067-0.0839]. 
With  the  above  values,  the  term  (24)  is  far  less  than  2LNQ,  which  allows  SINR  (23)  to  be 
very  close  to  its  upper  bound. 
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IV.  Conclusions 


In  this  report,  suppression  of  frequency  modulation  interferers  in  GPS  using  antenna 
arrays  and  subspace  projection  techniques  is  examined.  It  is  shown  that  both  the  spatial 
and  temporal  signatures  of  the  signals  impinging  on  the  multi-sensor  receiver  assume 
similar  roles  in  the  receiver  performance.  The  spatial  and  temporal  signatures’  respective 
correlation  coefficients  appear  as  square  multiplicative  products  in  the  SINR  expression. 
However,  due  to  the  Gold  code  structure  and  the  length  of  the  PN,  the  differences  in  the 
temporal  characteristics  of  the  jammer  and  the  C/A  code  yield  negligible  temporal 
correlation  coefficients.  These  small  values  allow  the  receiver  to  perform  very  close  to 
the  no-jamming  case,  irrespective  of  the  satellites  and  the  jammers'  angles  of  arrival.  The 
array  fundamental  offering  in  the  underlying  interference  suppression  problem  is  through 
its  gain,  which  is  determined  by  the  number  of  antennas  employed  at  the  receiver,  as 
shown  in  Fig.  1. 
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Chapter  4 


Array  Processing  for  Nonstationary  Interference 
Suppression  in  DS/SS  Communications  Using  Subspace 

Projection  Techniques 


Abstract 

Combined  spatial  and  time-frequency  signatures  of  signal  arrivals  at  a  multi-sensor  array 
are  used  for  nonstationary  interference  suppression  in  direct-sequence  spread-spectrum 
(DS/SS)  communications.  With  random  PN  spreading  code  and  deterministic 
nonstationary  interferers,  the  use  of  antenna  arrays  offers  increased  DS/SS  signal 
dimensionality  relative  to  the  interferers.  Interference  mitigation  through  spatio-temporal 
subspace  projection  technique  leads  to  reduced  DS/SS  signal  distortion  and  improved 
performance  over  the  case  of  a  single  antenna  receiver.  The  angular  separation  between 
the  interference  and  desired  signals  is  shown  to  play  a  fundamental  role  in  trading  off  the 
contribution  of  the  spatial  and  time-frequency  signatures  to  the  interference  mitigation 
process.  The  expressions  of  the  receiver  SINK  implementing  subspace  projections  are 
derived  and  numerical  results  are  provided. 
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1.  Introduction 


There  are  several  methods  that  have  been  proposed  for  interference  suppression  in 
DS/SS  communications,  most  have  been  related  to  one  domain  of  operation  [1],  [2].  These 
methods  include  the  narrowband  interference  waveform  estimation  [3],  frequency  domain 
interference  excision  [4],  zero-forcing  techniques  [5],  adaptive  subspace-based  techniques 
[6],  [7],  and  minimum-mean-square  error  (MMSE)  interference  mitigation  techniques  [8]. 

Nonstationary  interferers,  which  have  model  parameters  that  change  with  time,  are 
particularly  troublesome  due  to  the  inability  of  a  single  domain  mitigation  algorithm  to 
adequately  remove  their  effects.  The  recent  development  of  the  quadratic  time-frequency 
distributions  (TFDs)  for  improved  signal  power  localization  in  the  time-frequency  plane 
has  motivated  several  new  approaches  for  excision  of  interference  with  rapidly  time- varying 
frequency  characteristics  in  the  DS/SS  communication  systems.  Comprehensive  summary 
of  TFD-based  interference  excision  is  given  in  reference  [9].  The  two  basic  methods  for 
time-frequency  excision  are  based  on  notch  filtering  and  subspace  projections.  Utilization 
of  the  interference  instantaneous  frequency  (IF),  as  obtained  via  TFDs,  to  design  an  open 
loop  adaptive  notch  filter  in  the  temporal  domain,  has  been  thoroughly  discussed  in  [10], 

[11] .  Subspace  projection  methods,  commonly  used  for  mitigating  narrowband  interference 

[12] ,  [13],  have  been  recently  introduced  for  suppression  of  frequency  modulated  (FM) 
interference  and  shown  to  properly  handle  multi-component  interference,  reduce  the  self¬ 
noise,  and  improve  the  receiver  performance  beyond  that  offered  by  other  time-frequency 
based  techniques  [14],  [15],  [16]. 

The  main  purpose  of  this  chapter  is  to  integrate  spatial  and  temporal  processing  for 
suppression  of  nonstationary  interferers  in  DS/SS  systems.  Specifically,  we  extend  the 
projection-based  interference  mitigation  techniques  in  [14],  [15],  [16]  to  multi-sensor  array 
receivers.  The  proposed  multi-sensor  interference  excision  technique  builds  on  the  offerings 
of  quadratic  time- frequency  distributions  for  estimation  of  1)  the  time-frequency  subspace 
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and  time- frequency  signature  of  nonstationary  signals,  and  2)  the  spatial  signature  of 
nonstationary  sources  using  direction  finding  and  blind  source  separations.  With  the 
knowledge  of  the  time-frequency  and  spatial  signatures,  the  objective  is  to  effectively 
suppress  strong  nonstationary  interferers  with  few'  array  sensors.  The  proposed  technique 
does  not  require  the  knowledge  of  the  array  response  or  channel  estimation  of  the  DS/SS 
signal,  but  it  utilizes  the  distinction  in  both  of  its  spatial-  and  time-frequency  signatures 
from  those  of  the  interferers  that  impinge  on  the  array.  With  the  combined  spatial-time- 
frequency  signatures,  the  projection  of  the  data  vector  onto  the  subspace  orthogonal  to 
that  of  the  interferers  leads  to  improved  receiver  performance  over  that  obtained  using  the 
subspace  projection  in  the  single-sensor  case. 

The  rest  of  the  chapter  is  organized  as  follows.  In  Section  II,  the  signal  model  is 
described.  Section  III  briefly  reviews  the  subspace  projection  technique.  We  present 
in  Section  IV  blind  beamforming  based  on  subspace  projection  and  derive  the  receiver 
output  signal-to-interference-plus-noise  ratio  (SINK).  Several  numerical  results  are  given 
in  Section  V.  Section  VI  concludes  this  chapter. 

II.  Signal  Model 

In  DS/SS  communications,  each  symbol  is  spread  into  L  =  T/Tc  chips,  where  T  and 
Tc  are,  respectively,  the  symbol  duration  and  chip  duration.  We  use  discrete-time  form, 
where  all  signal  arrivals  are  sampled  at  the  chip-rate  of  the  DS/SS  signal.  The  symbol-rate 
source  signal  is  denoted  as  s(n),  and  the  aperiodic  binary  spreading  sequence  of  the  nth 
symbol  period  is  represented  by  c(n,  1)  G  ±1, 1  =  0, 1,  •  •  • ,  L  —  1.  The  chip-rate  sequence 
of  the  DS/SS  signal  can  be  expressed  as 

d{k)  —  s{n)c{n,  1)  with  k  —  nL  +  l.  (1) 

For  notation  simplicity,  we  use  c(/)  instead  of  c(n,  1)  for  the  spreading  sequence. 

We  consider  an  antenna  array  of  N  sensors.  The  propagation  delay  between  antenna 
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elements  is  assumed  to  be  small  relative  to  the  inverse  of  the  transmission  bandwidth,  so 
that  the  received  signal  at  the  .V  sensors  are  identical  to  within  complex  constants.  The 
received  signal  vector  of  the  DS/SS  signal  at  the  array  is  expressed  by  the  product  of  the 
chip-rate  sequence  d{k)  and  its  spatial  signature  h, 

x,(A-)  =  d{k)h.  (2) 

The  channel  is  restricted  to  flat-fading,  and  is  assumed  fixed  over  the  symbol  length,  and 
as  such  h  in  the  above  equation  is  not  a  function  of  k. 

The  array  vector  associated  with  a  total  of  V  interference  signals  is  given  by 

^uik)  =  (3) 

j=i 

where  a,  is  the  array  response  to  the  zth  interferer,  u,(A).  Without  loss  of  generality,  we 
set  ||h||^  =  N  and  ||a,||p  =  N,  i  =  l.2.  -  ,U,  where  ||  •  ||f  is  the  Frobenius  norm  of  a 
vector.  The  input  data  vector  is  the  sum  of  three  components, 

u 

x(A:)  =  x^ik)  +  Xu  (A)  +  h{k)  =  d{k)h  +  Y^a,Ui{k)  +  h{k)  (4) 

j=i 

where  b(A)  is  the  additive  noise  vector.  In  regards  to  the  above  equation,  we  make  the 
following  assumptions. 

Al)  The  information  symbols  s(n)  is  a  wide-sense  stationary  process  with  £’[s(n)s*(n)]  = 
1,  where  the  superscript  *  denotes  complex  conjugation.  The  spreading  sequence  c{k)  is  a 
binary  random  sequence  with  E[c{k)c{k  -I-  /)]  =  5{l),  where  6{l)  is  the  delta  function.^ 

A2)  The  noise  vector  b(A)  is  zero-mean,  temporally  and  spatially  white  with 

£’[b(A)b^(A  -|-  /)]  =  0  for  all  /, 

and 

E[h{k)h^  {k  +  1)]  =  a6{l)lN, 

'This  assumption  is  most  suitable  for  military  applications  and  P-code  GPS. 
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where  a  is  the  noise  power,  the  superscripts  ^  and  ^  denote  transpose  and  conjugate 
transpose,  respectively,  and  is  the  N  x  N  identity  matrix. 

A3)  The  signal  and  noise  are  statistically  uncorrelated. 

III.  Subspace  Projection 

The  aim  of  subspace  projection  techniques  is  to  remove  the  interference  components 
before  despreading  by  projecting  the  input  data  on  the  subspace  orthogonal  to  the  inter¬ 
ference  subspace,  as  illustrated  in  Fig.  1.  The  receiver  block  diagram  is  shown  in  Fig.  2. 

A  nonstationary  interference,  such  as  an  FM  signal,  often  shares  the  same  bandwidth 
with  the  DS/SS  signal  and  noise.  As  such,  for  a  chirp  signal  or  a  signal  with  high-order  fre¬ 
quency  laws,  the  signal  spectrum  may  span  the  entire  frequency  band,  and  the  sample  data 
matrix  loses  its  complex  exponential  structure  responsible  for  its  singularity.  Therefore, 
the  interference  subspace  can  no  longer  be  obtained  from  the  eigendecomposition  of  the 
sample  data  matrix  [12],  [14]  or  the  data  covariance  matrix  [13],  as  it  is  typically  the  case 
in  stationary  environments.  The  nonstationary  interference  subspace,  however,  may  be 
constructed  using  the  interference  time-frequency  signature.  Methods  for  estimating  the 
instantaneous  frequency,  instantaneous  bandwidth,  and  more  generally,  a  time-frequency 
subspace,  based  on  the  signal  time-frequency  localization  properties  are,  respectively,  dis¬ 
cussed  in  references  [17],  [18],  [14]. 

For  the  general  class  of  FM  signals,  and  providing  that  interference  suppression  is 
performed  separately  over  the  different  data  symbols,  the  interference  subspace  is  one¬ 
dimensional  in  an  L-dimensional  space.  We  note  that  since  an  FM  interference  has  a 
constant  amplitude,  its  respective  data  vector  can  be  determined  from  the  IF  up  to  a 
complex  multiplication  factor.  The  unit  norm  normalization  of  this  vector  represents  the 
one-dimensional  interference  subspace  basis  vector.  Among  candidate  methods  of  IF  es¬ 
timation  is  the  one  based  on  the  time-frequency  distributions.  For  example,  the  discrete 
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form  of  Cohen's  class  of  TFD  of  a  signal  3-(f)  is  given  by  [19] 

OC  OC 

Dxx{i-f)=  X]  X  <P{m.T)x{t  +  m  + T)x*{t  +  m  -  .  (5) 

m=—OQ  r— — OC 

where  (f){m,T)  is  a  time-frequency  kernel  that  could  be  signal-dependent.  The  TFD  con¬ 
centrates  the  interference  signal  pow’er  around  the  IF  and  makes  it  visible  in  the  noise 
and  PN  sequence  background  [17],  [20].  It  has  been  shown  that,  for  linear  FM  signals, 
Radon-Wigner  transform  provides  improved  IF  estimates  over  the  TFD  [21].  Parametric 
methods  using  autoregressive  model  have  also  been  proposed  [22]. 

Other  nonstationary  interference  with  instantaneous  bandwidth  or  spread  in  the  time- 
frequency  domain  are  captured  in  a  higher-dimension  subspace.  In  this  case,  the  inter¬ 
ference  subspace  can  be  constructed  from  the  interference  localization  region  fl  in  the 
time-frequency  domain  (see,  for  example,  [14]).  The  subspace  of  interest  becomes  that 
which  fills  out  the  interference  time-frequency  region  D  energetically,  but  has  little  or  no 
energy  outside  f2. 

Interference-free  DS/SS  signals  are  obtained  by  projecting  the  received  data  vector  (in 
the  temporal  domain  processing,  the  vector  consists  of  data  samples  at  different  snapshots) 
on  the  subspace  orthogonal  to  the  interference  subspace. 

A.  Temporal  Processing 

In  the  single-sensor  receiver,  the  input  data  is  expressed  as 

u 

x{k)  =  Xs{k)  +  Xuik)  +  b{k)  =  d{k)  +  X^j(^)  +  (^) 

7  =  1 

Using  L  sequential  chip-rate  samples  of  one  symbol  of  the  received  signals  at  time  index 

k,  we  obtain  the  following  input  vector 

[x{k)  x{k  —  1)  ■■■  x{k  —  L-{-l)]^ 

^[x,{k)  Xs{k-1)  •••  Xs{k-LAl)f 
+  [xu{k)  Xu{k-1)  ■  ■  ■  Xu{k  -  L  +  1)]^ 

+  [b{k)  b{k-l)  ■■■  b{k-L  +  l)f 
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or  simply 


A'(/c)  =  Xs{k)  +  X^ik)  +  B{k).  (8) 

We  drop  the  variable  k  for  simplicity,  with  the  understanding  that  processing  is  per¬ 
formed  over  the  nth  symbol  that  starts  at  the  A:th  chip.  Then,  equation  (8)  becomes 

X  =  Xs  +  Xu  +  B.  (9) 

Below,  we  relax  the  FM  condition  used  in  [12],  [15]  that  translates  to  a  single  dimension 
interference.  The  general  case  of  an  interference  occupying  higher  dimension  subspace  is 
considered.  We  assume  that  the  ith  interferer  spans  A/,  dimensional  subspace,  defined  by 
the  orthonormal  basis  vectors,  1^,1,  Vi.2,  •  •  • ,  and  the  different  interference  subspaces 
are  disjoint.  Define 

l^  =  [niH2  •••  Vod  (10) 

and  let  M  =  A/,  as  the  number  of  total  dimensions  of  the  interferers.  W’ith  L  >  A/, 

the  L  X  M  matrix 


V  =  [\\V■^  ■■■Vu],  =  ^  (11) 

is  full  rank  and  its  columns  span  the  combined  interference  subspace  J.  The  respective 
projection  matrix  is 

P  =  v(y^VyW^.  (12) 

The  projection  matrix  associated  with  the  interference  orthogonal  subspace,  G,  is  then 
given  by  ^ 

p  =  1l-  v{v^vyH^^.  (13) 

When  applied  to  A”,  matrix  P  projects  the  input  data  vector  onto  G,  and  results  in 

Xi.  =  pa:  =  PXs  +  PB,  (14) 

which  no  longer  includes  any  interference  component. 
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The  single-sensor  receiver  implementing  siibspace  projection  for  excision  of  a  single 
instantaneously  narrowband  FM  interferer  (i.e..  U  =  1.  iMi  =  1)  in  DS/SS  communications 
is  derived  in  [23].  The  receiver  SINK  is  shown  to  be 


SINK  = 


{L-1? 


L-  1 


For  typical  values  of  L,  {L  -  2)/{L  -  1)  ^  1,  and  equation  (15)  can  be  simplified  as 


SINK  « 


L  -  1 
a  +  1/L 


(16) 


Compared  to  the  interference-free  environment,  where  the  receiver  SINK  is  Lja,  nonsta¬ 
tionary  interference  suppression  in  (16)  is  achieved  by  reducing  the  processing  gain  by  1 
and  increasing  the  noise  power  by  the  self-noise  factor  of  1/L. 


IV.  Subspace  Projection  in  Multi-Sensor  Receiver 


In  this  section,  we  consider  nonstationary  interference  excision  in  multi-sensor  receivers 
using  subspace  projections.  We  note  that  if  the  subspace  projection  method  discussed  in 
Section  III  is  extended  to  an  V-element  array  by  suppressing  the  interference  indepen¬ 
dently  in  each  sensor  data  and  then  combining  the  results  by  maximum  ratio  combining 
(see  Fig.  3),  then  it  is  straightforward  to  show  that  the  receiver  SINK  is  given  by 


SINK 


N[L-\) 
a  -h  N/L 


(17) 


The  above  extension,  although  clearly  improves  over  (16),  does  not  utilize  the  potential 
difference  in  the  spatial  signatures  of  signal  arrivals,  and,  therefore,  is  inferior  to  the 
receiver  proposed  in  this  Section. 


A.  Spatio-Temporal  Signal  Subspace  Estimation 

To  construct  the  spatio-temporal  signal  subspace  of  the  interference  signals,  it  is  impor¬ 
tant  to  estimate  both  the  time-frequency  signature  (or  subspace)  and  the  spatial  signature 
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of  each  interferer.  The  IF  estimation  of  an  FM  interference  signal  based  on  time-frequency 
distribution  is  addressed  in  Section  III.  It  is  noteworthy  that  when  multiple  antennas  are 
available,  the  TFD  may  be  computed  at  each  sensor  data  separately  and  then  averaged 
over  the  array.  This  method  has  been  shown  in  [24]  to  improve  the  IF  estimation,  as  it 
reduces  noise  and  crossterms  that  often  obscure  the  source  true  power  localization  in  the 
time- frequency  domain. 

On  the  other  hand,  the  estimation  of  source  spatial  signature  can  be  achieved,  for  ex¬ 
ample,  by  using  direction  finding  and  source  separation  techniques.  When  the  interference 
signals  have  clear  bearings,  methods  like  MUSIC  [25]  and  maximum  likelihood  (ML)  [26] 
can  be  used  to  estimate  the  steering  matrix  of  the  interference  signals.  These  methods 
can  be  revised  to  incorporate  the  TFD  of  the  signal  arrivals  for  improved  performance 
[27],  [28].  On  the  other  hand,  in  fading  channels  where  the  steering  vector  loses  its  known 
structure  due  to  multipath,  blind  source  separation  methods  should  be  used  [29],  [30], 
[31].  Since  the  interferers  in  DS/SS  communications  often  have  relatively  high  power, 
good  spatial  signature  estimation  is  expected. 

More  conveniently,  the  spatial  signatures  can  be  simply  estimated  by  using  matched 
filtering  once  the  time- frequency  signatures  are  provided.  The  maximum  likelihood  esti¬ 
mator  for  the  vector  a^  is  obtained  as 

ii  =  H  (^)x(^)/i  I  H  u* {k)x{k)  I \f,  (18) 

k=0  fc=0 

where  Ui{k)  is  the  estimated  waveform  of  the  zth  interferer.  It  is  noted  that  the  possible 
phase  ambiguity  in  the  waveform  estimation  of  Ui{k)  does  not  affect  the  estimation  of  the 
spatial  signature.  For  slowly  varying  channels,  the  above  average  can  also  be  performed 
over  multiple  symbols  to  improve  the  estimation  accuracy. 

In  the  analysis  presented  herein,  we  assume  knowledge  of  the  interference  subspace  and 
its  angle-of-arrival  (AOA)  to  derive  the  receiver  SINR. 
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D.  Proposed  Technique 


The  subspace  projection  problem  for  nonstationary  interference  suppression  in  DS/SS 
communications  is  now  considered  within  the  context  of  multi-sensor  array  using  N  array 
elements.  We  use  one  symbol  DS/SS  signal  duration  (i.e.,  L  chip-rate  temporal  snapshots), 
and  stack  L  discrete  observations  to  construct  an  NL  x  1  vector  of  the  received  signal 
sequence  in  the  joint  spatio-temporal  domain.  In  this  case,  the  received  signal  vector  in 
(4)  becomes 


(19) 


or  simply 


(k)  —  1) 

... 

{k  —  L  1) 

= 

[xf  (k) 

^J{k 

-1)  •• 

■  xKA:- 

-i  +  i)]’ 

+  1 

xI(A-) 

Kl{k 

-1)  ••• 

^lik- 

i+Df 

+  1 

Vik) 

h^{k 

-1)  ... 

h^{k  - 

i  +  Df 

X  = 

X,-f  x„ 

+  B, 

(20) 


where  again  the  variable  h  is  dropped  for  simplicity. 

In  (19),  the  interference  vector  in  the  single-sensor  problem,  given  by  (7),  is  extended 
to  a  higher  dimension.  With  the  inclusion  of  both  temporal  and  spatial  samples,  the  mth 
basis  of  the  zth  interference  becomes 


—  U.m  ®  §i 

(21) 

V.  =  |V,,,  V,,,  . . .  V,„,] . 

(22) 

where  ®  denotes  the  Kronecker  product.  The  columns  of  the  LN  x  M  matrix 

V-[Vi  V2  •••  Vj,-]  (23) 


spans  the  overall  interference  signal  subspace. 

For  independent  spatial  signatures,  the  matrix  rank  is  M.  The  orthogonal  projection 
matrix  is  given  by 

P  =  I^^^-V(V"V)“V^.  (24) 
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The  projection  of  the  signal  vector  on  the  orthogonal  subspace  of  the  interferers’  yields 

Xi  =PX  =  PX, +  PB.  (25) 

The  block  diagram  of  the  proposed  method  is  presented  in  Fig.  4.  As  shown  in  the 
next  section,  effective  interference  suppression  can  be  achieved  solely  based  on  the  spatial 
signatures  or  the  time-frequency  signatures,  or  it  may  require  both  information. 

C.  Performance  Analysis 

Below  we  consider  the  performance  of  the  multi-sensor  receiver  system  implementing 
subspace  projections.  Recall  that 

=  0  for  any  i,  m  #  j,  n.  (26) 

and 

=  NIm,  (27) 

the  projection  matrix  P  becomes 

(28) 

The  signal  vector  X^  can  be  rewritten  as 

X,  =  [xj {k)  xj{k-l)  xj{k-  L  +  1)] 

=  [d{k)h^  d{k  -  l)h^  •  •  •  d{k-L  +  l)h^]^ 

=  s{n)  [c{L  -  1)  c(L  -  2)  •  •  •  c(0)]^  ®  h 
A  s(n)q, 

where  the  LN  x  1  vector 

q  =  [c(L  -  1)  c(L  -  2)  •••  c(0)f  ®h  Ac0h  (30) 

defines  the  spatio-temporal  signature  of  the  desired  DS/SS  signal,  q  is  the  extension  of 
the  DS/SS  code  by  replicating  it  with  weights  defined  by  the  signal  spatial  signature. 
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By  performing  despreading  and  beamforming,  the  symbol-rate  decision  variable  is  given 

b}' 


y{n)  =  q"X_L(A:)  =  s(n)q"Pq -f  q^PB(/c)  A  yi(n)  -f  y2{n) 


(31) 


where  yi{n)  is  the  contribution  of  the  desired  DS/SS  signal  to  the  decision  variable,  and 
y2{n)  is  the  respective  contribution  from  the  noise. 

The  SINK  of  the  array  output  becomes  (see  Appendix  A) 

2 

E^[y(k)] 


SINK  = 


(L-'£Am\ 
\  1=1 


var[y{k)] 


u_ 

>_ 

1=1 


Ej\/,is,p)  -  2y;{.iar' )  +  ^  (i- 

i=l 


N 


i=l 


(32) 


where  is  defined  in  (A. 9),  and  is  the  spatial  correlation  coefficient  between  the  spatial 
signatures  h  and  a,,  f  =  1,  2,  •  •  • .  [/,  and  is  given  by 

1 


(33) 


Note  that  when  the  noise  power  is  small,  i.e.,  a  1.  the  variance  of  yi  becomes  dominant, 
and  the  output  SINK  reaches  the  following  upper  bound 

u 


SINK 


high  SNR 


\  j=i 


(34) 


>^1=1  /  i=l 

This  result  is  affected  by  the  factors  L,  A/,,  |/?,|,  and  ^i.i  =  I,  -  ■  ■  ,U.  On  the  other  hand, 
when  the  noise  level  is  very  high,  i.e.,  a  ^  1,  the  noise  variance  plays  a  key  role  in 

determining  var[y{k)],  and  the  output  SINK  becomes 

2 


SINK 


low  SNR 


V  ^  b-EM.ia 


a 

A 


(35) 


i=l 
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Unlike  the  high  input  SNR  case,  the  output  SINK  in  (35)  also  depends  on  both  K  and  a. 
Comparing  (34)  and  (35).  it  is  clear  that  the  improvement  in  the  receiver  SINK  becomes 
more  significant  when  the  spatial  signatures  produce  small  spatial  correlation  coefficients 
and  under  high  SNR. 

Next,  we  consider  some  specific  important  cases.  When  Pi  —  0,i  =  1,  -  ■  ■  ,U,  t'ar[yi(n)]  = 
0,  the  receiver  SINR  in  (32)  becomes  SINR  =  LN/a.  This  is  to  say,  the  output  SINR 
is  improved  by  a  factor  of  LA’’  over  the  input  signal-to-noise  ratio  (SNR)  (not  the  input 
SINR!).  This  implies  that  the  interferers  are  suppressed  by  spatial  selectivity  of  the  array 
and  their  suppression  does  not  cause  any  distortion  of  the  temporal  characteristics  of  the 
DS/SS  signal.  The  DS/SS  signal  in  this  case  enjoys  the  array  gain  that  contributes  the 
factor  N  to  the  SINR. 


For  a  single  FM  interferer  {U  =  1,  Mi  =  1),  equation  (32)  becomes 

SINR  =  ^ - . 


(N _ 

1 

2 

1 

(-1) 

h-iAih 

It  is  easy  to  show  that  SINR  in  (36)  monotonously  decreases  as  |/?i|  increases,  and  the 
lower  bound  of  the  SINR  is  reached  for  Pi  =  1,  which  is  the  case  of  the  desired  DS/SS 
signal  and  the  interference  signal  arriving  from  the  same  direction.  With  a  unit  value  of 

SINR  =  7 - (37) 


/  2\  cr  N  .  ' 

(l-E)  +  jf(i-l)  T  +  - 

This  result  is  the  same  as  that  of  the  single-sensor  case  developed  in  [15],  except  for  the 
appearance  of  the  array  gain,  N,  for  the  desired  DS/SS  signal  over  the  noise.  This  equa¬ 
tion  also  coincides  with  (17).  That  is,  the  independent  multi-sensor  subspace  projection, 
illustrated  in  Fig.  3,  results  in  the  same  output  SINR  with  the  proposed  multi-sensor 
subspace  projection  method  when  |/?il  =  0. 

On  the  other  hand,  the  maximum  value  in  (36)  corresponds  to  P  —  0,  and  is  equal  to 
SINR  =  LN/a,  as  discussed  above.  For  the  illustration  of  the  SINR  behaviour,  we  plot 
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in  Fig.  5  the  SINK  in  (36)  versus  |/?]|  for  a  two-sensor  array,  where  L  =  64,  and  one  FM 
jammer  is  considered  with  M  —  7.  The  input  SXR  is  OdB. 

Given  the  temporal  and  spatial  signatures,  the  proposed  technique  simplifies  to  two 
consecutive  tasks.  The  first  is  to  estimate  the  spatio-temporal  signature.  When  using 
multiple  antenna  receivers,  a  basis  vector  of  the  orthogonal  projection  matrix  is  obtained 
by  the  Kronecker  product  of  a  jammer’s  temporal  signature  and  its  spatial  signature,  that 
results  in  the  LN  x  LN  orthogonal  project  matrix  instead  of  L  x  L  in  the  single  antenna 
case.  The  second  task  is  jammer  suppression  via  subspace  projection.  This  involves  the 
multiplication  of  an  LN  x  LN  matrix  and  an  LN  x  1  vector. 

Note  such  increase  in  computations  is  natural  due  to  increase  of  dimensionality.  It 
is  noteworthy  that  array  processing  expands  overall  space  dimensionality  but  maintains 
the  jammer  subspace  dimension.  As  a  result,  it  yields  improved  SINR  performance  over 
temporal  processing  or  spatial  processing  only  methods. 

Y.  Numerical  Results 

A  two-element  array  is  considered  with  half-wavelength  spacing.  The  DS/SS  signal  uses 
random  spreading  sequence  with  L  —  64.  The  AOA  of  the  DS/SS  signal  is  0  degree  from 
broadside  {6d  =  0°). 

We  consider  twm  interference  signals.  Each  interference  signal  is  assumed  to  be  made  up 
of  uncorrelated  FM  component  with  Mj  —  7,  z  =  1,2.  The  overall  interference  subspace 
is  M-14.  The  AOAs  of  the  two  interferers  are  Oj  =  [40°,  60°].  The  respective  spatial 
correlations  in  this  example  are  |^i|  =  0.53  and  |/?2|  =  0.21.  Note  that,  in  the  subspace 
projection  method,  the  output  SINR  is  independent  of  the  input  jammer-to-signal  ratio 
(JSR),  since  the  interferers  are  entirely  suppressed,  regardless  of  their  power.  Fig.  6  shows 
the  receiver  SINR  versus  the  input  SNR.  The  upper  bounds  correspond  to  interference-free 
data.  For  high  input  SNR,  the  receiver  SINR  is  decided  by  the  induced  signal  distortion, 
described  by  the  variance  given  in  (A.  10).  It  is  evident  from  Fig.  6  that  the  two-antenna 
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receiver  outperforms  the  single-antenna  receiver  case  by  a  factor  much  larger  than  the 
array  gain.  Since  the  output  SINK  in  the  two- antenna  receiver  highly  depends  on  the 
spatial  correlation  coefficients,  the  curves  corresponding  to  a  two-sensor  array  in  Fig.  6 
will  assume  different  values  upon  changing  or/and  (32-  The  best  performance  is  achieved 
at  —  132  =  0. 

Fig.  7  shows  the  receiver  SINK  versus  the  number  of  chips  per  symbol  (L).  We  let  L 
vary  from  8  to  4096,  whereas  the  input  SNR  is  fixed  at  0  dB.  The  two  interference  signals 
are  incident  on  the  array  with  angles  6j  -  [40°,  60°].  They  are  assumed  to  maintain  their 
time-frequency  spread  with  increased  value  of  L.  As  such,  the  respective  dimentions  of 
their  subspaces  grow  propotional  to  the  number  of  chips  per  symbol.  In  this  example,  the 
dimension  of  each  interference  signal  is  assumed  to  be  10  percent  of  L  (round  to  the  nearest 
integer).  The  ouput  SINK  improvement  by  performing  array  processing  at  different  L  is 
evident  from  this  figure.  It  is  seen  that,  unlike  the  case  of  the  instantaneously  narrowband 
FM  interference,  where  the  output  SINK  increases  rapidly  as  L  increases,  the  output  SINK 
in  the  underlying  scenario  ceases  to  increase  as  L  assumes  large  values.  This  is  because 
the  rank  of  the  interference  signal  subspace  increases  with  L. 

In  Fig.  8  we  investigate  the  receiver  SINK  performance  versus  the  number  of  array 
sensors.  In  this  figure,  L  is  set  at  64,  and  the  input  SNR  is  0  dB.  Two  interference  sig¬ 
nals  composed  of  uncorrelated  FM  components  are  considered,  and  Mj  =  7,  f  =  1,2,  are 
assumed.  Two  examples  are  used  to  examine  the  effect  of  different  AOAs.  In  the  first  ex¬ 
ample,  Oj  -  [40°,  60°].  The  output  SINR  improves  sharply  as  the  number  of  array  sensors 
increases  from  one  to  three,  beyond  which  the  improvement  becomes  insignificant.  The 
differences  in  the  above  AOAs  of  the  desired  DS/SS  signal  and  the  interference  signals 
are  relatively  large,  and  a  small  number  of  array  sensors  leads  to  negligible  spatial  cor¬ 
relation  coefficients.  We  also  show  a  case  with  closely  spaced  interference  signals  where 
Oj  =  [5°,  15°].  In  this  case,  the  output  SINR  slowly  improves  as  the  number  of  array 
sensors  increases. 
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It  is  noted  that,  when  we  consider  a  specific  case,  the  output  SINK  does  not  increase 
monotonously  with  the  number  of  array  sensors.  This  is  because  the  relationship  between 
the  spatial  correlation  coefficient  and  the  AOAs  is  by  itself  not  monotonous.  Nevertheless, 
when  we  consider  the  general  case  with  different  AOA  combinations,  high  number  of  array 
sensors  often  reduce  the  spatial  correlation  coefficients. 

VI.  Conclusions 

In  this  chapter,  subspace  projection  techniques  were  employed  to  suppress  nonstationary 
interferers  in  direct  sequence  spread  spectrum  (DS/SS)  communication  systems.  Interfer¬ 
ence  suppression  is  based  on  the  knowledge  of  both  the  interference  time-frequency  and 
spatial  signatures.  While  the  former  is  based  on  instantaneous  frequency  information  that 
can  be  gained  using  several  methods,  including  time-frequency  distributions,  the  later  can 
be  provided  from  applying  higher  resolution  methods  or  blind  source  separation  techniques 
to  the  signal  arrivals. 

The  differences  between  the  DS/SS  signal  and  interference  signatures  both  in  the  time- 
frequency  and  spatial  domains  equip  the  projection  techniques  with  the  ability  to  remove 
the  interference  with  a  minimum  distortion  of  the  desired  signal. 

The  receiver  performance  based  on  subspace  projections  was  analyzed.  It  was  shown 
that  the  lower  performance  bound  is  obtained  when  the  sources  have  the  same  angular 
position.  In  this  case,  the  problem  becomes  equivalent  to  a  single-antenna  receiver  with 
only  the  presence  of  the  array  gain.  On  the  other  hand,  the  upper  bound  on  performance 
is  reached  in  the  interference-free  environment  and  also  corresponds  to  the  case  in  which 
the  spatial  signature  of  the  interference  is  orthogonal  to  that  of  the  DS/SS  signal. 

Numerical  results  were  presented  to  illustrate  the  receiver  SINK  dependency  on  spatial 
correlation  coefficient,  input  SNR,  and  the  PN  sequence  length. 
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Appendix  A 


To  derive  the  output  SINK  expression,  we  use  s{n)  =  +1  (the  output  SINK  is  indepen¬ 
dent  of  s{n)  and  same  result  follows  when  s(n)  =  -1).  Then, 

£|!/.(»)1  =£[q"Pq] 

=  £[q''(l-ivv«)q] 

=  £[q''q]-;i£[q''W»ql 

=  LAT  -  q"  I:  V,,  f:  V"q 

L  n=i  <2=1 

1  r  U  Mi  U  Mi  1 


1  U  iVit  U 

=  LN  -  iiS  q"  E  Y.  Vi,.™.  E  E  V,"„.q  ■ 


ii=lmi=l  12  =  1  mo  =  l 


It  is  straightforward  to  show  that 


q''Vi„  =  (c  ®  h)"  {v;,„  ®  a.)  =  (s’'V',„)  ®  (h''a,)  =  A'ft  Y  VUW)-  (A.2) 


Using  the  orthogonal  property  of  the  spreading  sequence  Al),  (A.l)  becomes 

V  ^'^>1  L-l  V  ^^'2  L-l 

E[y,{n)]  ^LN-NE  E  E  (^iM/i)  E  E  E  K;m2 ^2) 

ii=l  mi=l/i=0  12=1  7712=1/2=0 


7711  =  1  /i=0 

LN  -  iv|:  ^  |V^m(0Pc2(/) 

t=l  m—l  1=0 


\  i=l 


(A.3) 


Due  to  the  zero-mean  property  of  noise  (assumption  A2),  E[y2{n)\  =  0.  Accordingly, 


f;[y(n)]  =  iE[yi(n)]  =  iV 


(A.4) 


It  is  clear  from  (A.4)  that  the  increase  in  the  space  dimensionality  from  L  to  NL  does  not 
simply  translate  into  a  corresponding  increase  in  the  desired  mean  value,  or  subsequently 
in  the  processing  gain.  Also,  from  assumption  A3),  the  cross-correlation  between  yi(n) 
and  y2(n)  is  zero,  i.e., 

£:[yi(n)y2(n)]  =  E[yi{n)yl{n)\  =  0.  (A.5) 
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Therefore,  the  mean  square  value  of  the  decision  variable  is  made  up  of  only  two  terms. 


E  [ly(tOi^]  =  E  [|yi(n)|^]  +  E  [\y2{n)\^  . 


(A.6) 


The  first  term  is  the  mean  square  value  of  yi(n).  From  (26),  we  have 
E  [|yi(n)i^]  =  E  [q^Pqq^P^q] 
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In  practice,  7,;  takes  negligible  values,  and  equation  (A. 7)  can  be  simplified  to 


(A.7) 


(A.8) 


(A.9) 


E  [|yi(n)|"]  =  {LNf  -  2LA^^  ^  M,|Ap  +  A^^  j  -  2E6iAI'j  •  (A.IO) 

The  value  of  depends  on  the  type  of  interference  signals.  Specifically,  when  the  ith 
interference  signal  is  made  up  of  a  single  FM  or  a  number  of  uncorrelated  FM  signal 
components,  then  the  basis  vectors  are  of  constant  modulus,  and 


c  -Mi 
L  ■ 

The  second  term  of  (A. 6)  is  the  mean-square  value  of  y2(tt), 

£:[!y2(n)|']  =.F[q^PB(/c)B"(fc)P^q] 

-  aE  [q^PP^q]  =  aE  [q"Pq]  =  aN  . 

The  variance  of  y(n)  is  given  by 
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Equation  (32)  follows  by  using  the  results  of  (A.4)  and  (A.  13). 


(A.13) 
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Fig.  1  Jammer  suppression  by  subspace  projection. 


Fig.  2  Block  diagram  of  single-sensor  subspace  projection. 
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Fig.  3  Block  diagram  of  independent  multi-sensor  subspace  projection. 


Fig.  4  Block  diagram  of  proposed  multi-sensor  subspace  projection. 
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